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Tidal Streams as Probes of the Galactic Potential 

Kathryn V. Johnston 1 , HongSheng Zhao 2 , David N. Spergel 3 and Lars Hernquist 4 

ABSTRACT 

We explore the use of tidal streams from Galactic satellites to recover the 
potential of the Milky Way. Our study is motivated both by the discovery of 



the first lengthy stellar stream in the halo ( |lrwin fc Totten 1998| ) and by the 



Q ! prospect of measuring proper motions of stars brighter than 20th magnitude in 



such a stream with an accuracy of ~ 4/zas/yr, as will be possible with the Space 
Interferometry Mission (SIM). We assume that the heliocentric radial velocities 
! of these stars can be determined from supporting ground-based spectroscopic 

surveys, and that the mass and phase-space coordinates of the Galactic satellite 
t^J- ! with which they are associated will also be known to SIM accuracy. Using 

results from numerical simulations as trial data sets, we find that, if we assume 
the correct form for the Galactic potential, we can predict the distances to the 
! stars as a consequence of the narrow distribution of energy expected along the 

r-| ! streams. We develop an algorithm to evaluate the accuracy of any adopted 

potential by requiring that the satellite and stars recombine within a Galactic 
lifetime when their current phase-space coordinates are integrated backwards. 
When applied to a four- dimensional grid of triaxial logarithmic potentials, with 
varying circular velocities, axis ratios and orientation of the major-axis in the 
disk plane, the algorithm can recover the parameters used for the Milky Way 
in a simulated data set to within a few percent using only 100 stars in a tidal 
stream. 

Subject headings: Galaxy: fundamental parameters, halo, kinematics 
and dynamics, structure 
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1. Introduction 

Tidal streams in the Galactic halo are a natural by-product of hierarchical 
structure formation, in which galaxies build up their mass by accreting less 
massive satellites. These features can be produced when matter, i.e. stars 
and/or gas, is liberated from a companion either by tidal shocking at the 
pericenter of the satellite's orbit or, in the case of globular clusters, through the 
evaporation of stars across the tidal boundary imposed by the Milky Way. The 
stripped material populates leading and trailing tidal streams that align with 
the orbit of the satellite for timescales comparable to or greater than the age 
of the Galaxy (e.g. ph, Lin fc Aarseth 1995| ; Grillmair et al. 1995| ; ]Iohnston"^ 



Spergel fc Hernquist 19951 ; Johnston, Hernquist & Bolte 1996, hereafter |JHB|) . 

The notion of using tidal streams as Galactic potentiometers has been 
considered previously by a number of authors ( [Lynden-Bell 1982] ; [Kuhn 1993| 



Grillmair 1998| ; |Zhao 1998| ). By assuming that several of the dwarf spheroidal 



satellites are tidal debris, Lynden-Bell (1982) was able to obtain an estimate of 
the mass of the Milky Way. Similarly, if the Magellanic Stream consists of gas 
tidally stripped from the Large and Small Magellanic Clouds, it can be used to 
constrain the potential flMurai fc Fujimoto 1980| ; [Lin fc Lynden-Bell 1982| ; |Lin 



Jones fc Klemola 1995|) . While these results are not definitive measurements 



owing to the controversial nature of the Magellanic Stream ( [Moore fc Davie 
1994|) , they demonstrate the power of this approach. 



There is now growing observational evidence for the existence of stellar 
tidal streams in the halo. Several globular clusters are known to possess 
excess unbound stars outside their tidal radii ( |Grillmair et al. 1995| |lrwin fc 



Hatzidimitriou 1995|; [Kuhn, Smith fc Hawley 1996|) and there are also moving 



groups in the halo with no obvious bound counterparts (for a review see 
Majewski, Hawley fc Munn 1996|) . On larger scales, the debris associated with 



the dwarf galaxy recently discovered in the constellation Sagittarius ( |Ibata,| 
Gilmorc fc Irwin 1994), has now been identified in horizontal branch (HB) and 



giant branch (GB) stars over 22° along a Galactic great circle roughly coinciding 
with the I = 0° plane ( Mateo et al. 1996; ; |Alard 1996| ; [Fahlman et al. 1996| ; 



Ibata et al. 1997] ). Moreover, Irwin & Totten's (1998) discovery of a carbon star 



trail encircling the Galaxy and aligned with the same plane provides the first 
example of a data set that samples the entire length of a stellar tidal stream. 

Upcoming satellite missions promise to provide accurate phase-space 
coordinates for individual stars in tidal streams. NASA's Space Interferometric 
Mission (SIM), scheduled for 2006, is a pointed instrument that will detect 
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stars as faint as 20th magnitude with accuracies of a few fias. ESA's Global 
Astrometric Interferometer for Astrophysics (GAIA) will survey more than a 
billion stars across the entire sky with an astrometric precision of < 10/ias. 
This represents an improvement over results obtained with the HIPPARCOS 
satellite by a factor of about a thousand in accuracy and more than a million 
in the volume sampled, and makes the idea of measuring the proper motions 
of individual HB and GB stars (and hence carbon stars) at distances < 100 
kpc feasible. At distances of tens of kiloparsecs, these future missions, when 
supplemented with radial velocities, will accurately measure five out of the six 
phase-space coordinates of a star. 

In this Letter, we investigate the extent to which the potential of the Galaxy 
can be recovered using a data set such as Irwin & Totten's (1998) carbon star 
stream, assuming that phase-space positions can be inferred with the accuracy 
of the SIM satellite. At the appropriate distances, parallaxes measured by 
SIM (~ Afias) will be the least well-determined phase-space component. The 
implied uncertainty in distances, ~ (-D/20kpc) 2 kpc for stars at a distance D, 
will exceed the thickness of the tidal stream. In §3.1, we show that using energy 
conservation in the correct background potential, together with the other five 
phase-space coordinates, yields a more precise distance estimate. We describe 
simulations that provide trial data sets for our proposed methods in §2, and use 
them to estimate the accuracy of distances to debris stars constrained by energy 
conservation in §3.1. We apply our algorithm for recovering the potential of the 
Milky Way in §3.2, and summarize conclusions in §4. 



2. Trial Data Sets and Galactic Model 

This paper is based on an analysis of simulated debris trails produced by 
tidal disruption, as reported in |JHB| . In these simulations, the Milky Way 
potential $mw was taken to be smooth, static and axisymmetric, and was 
represented by a three-component bulge-disk-halo model. In each simulation, 
the satellite was represented by 10 4 particles, whose self-gravity was calculated 
using a self-consistent field code ( |Hernquist fc Ostriker 199^ ). The particles 



were initially distributed as a Plummer model, and evolved in the Milky Way 



potential for 10 Gyrs. The reader is referred to |JHB| for full details of the 
simulations. 

We tested our algorithm on all the models (1-12) in JHB and, in the figures 
which follow, we illustrate the results using their Model 11. This satellite had 
an initial mass 2.893 x 1O 7 M and scale-length 0.602kpc. It was on an orbit 
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inclined at 50° relative to the Galactic disk with an orbital period ~ 2Gyrs and 
peri- and apo-center distances of ~ 40kpc and ~ 160kpc. We also performed 
two new simulations with the same initial conditions as in Model 11, but with 
the spherical halo component of the Milky Way potential replaced by oblate and 
triaxial components of the form 

$haio(x, y,z) = V -f^ ln(x 2 + y 2 /p 2 + z 2 /q 2 + c 2 ). (1) 

In the remainder of the paper we shall refer to these new simulations as Models 

13 and 14. In all cases v C i rc = 181km/s and c = 12kpc. In Models 1-12 

0, q) = (1, 1), in Model 13 (p, q) = (1, 0.75) and in Model 14 (p, q) = (0.95, 0.75). 



3. Results 

3.1. Energy Distances 

Tidal debris tends to become unbound from a satellite of mass m sat located 
at distance R from the center of the Galaxy on the physical scale given by the 
tidal radius 

'- = *(^)* (2) 

( [King 1962| ), where Mr is the mass of the Galaxy enclosed within R. Hence the 
orbital energies E of material in the debris will lie approximately within 

d$ _ GM R 

e — r tide "TT^ — '"tide — 

dR R A 



of the satellite's orbital energy, E sat QTremaine 1993| , |Johnston 1998| ). Equations 



(0) and (H) should be evaluated not at the current position of the satellite, 
but at the pericenter of its orbit since most of the mass loss will occur where 
the tidal field of the Milky Way is strongest. The top panel of Figure |I] shows 
E for all particles which were no longer bound to the satellite at the end of 
the simulation of Model 11, with the left hand axis in physical units and the 
right hand axis in units relative to E 8at and scaled by e given in equation (H). 
Note that the debris divides into two populations, corresponding to stars in the 
leading/trailing stream having E offset by =F5e/4 from E sat . Each population 
has a width ~ 3e/4 about these average offsets. 

Even with the astrometric accuracy of SIM, we will be able to measure 
distances only to a precision ^10% beyond 20kpc. However, if we can identify a 
population of stars which were once associated with a satellite and we assume 
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a form $mw f° r the Galactic potential, it is possible to estimate their distances 
from the expected energy distribution in the satellite's debris. To test this idea 
we "observe" the angular position (l,b), line of sight velocity t>i os and proper 
motion (/^,/^,) of each of the unbound particles at the end of our simulations 
from an assumed heliocentric viewpoint in the disk plane 8.5 kpc from the center 
of the Galaxy, where the motions are defined in a Galactic rest frame. Each 
particle is assigned E = E slA =p 5e/4 if it is ahead/behind the satellite along the 
orbit, and its heliocentric distance d is then calculated by solving the equation 

E = \ [vL + d 2 (tf + fil)] + $ MW (d, I, b). (4) 

The bottom panel of Figure [I] plots the error in the distance estimated for 
particles in Model 11 using this method. As can be seen from the right hand 
axis of this figure, we expect this distance estimate to be good to within a few 
r t idc- In contrast, the solid lines show the much lower accuracy of distances 
measured from parallaxes with the SIM satellite (as described in §1). 



3.2. Constraining the Galactic Potential 

If we observe a set of stars that were once associated with a satellite of 
known mass and phase-space coordinates, then they must be on orbits whose 
trajectories intersect with the satellite within the lifetime of the Galaxy. The 
following algorithm assigns a quality factor for each assumed form $>(x,y, z) for 
the Galactic potential based on how well it can extrapolate the orbits of tidal 
debris stars back into the progenitor satellite galaxy. 

1. For each assumed potential we integrate the satellite's orbit backwards and 
calculate r ti d e (eqn 0) and the energy scale e (eqn ||) at peri center. 

2. For each star in the debris with measured proper motion, angular 
position and radial velocity we: (i) create n tcs t particles with energies 
E evenly distributed in the range ±3e/4 about (E sat =p 5e/4) if the 
star is ahead/behind the satellite along its orbit; (ii) estimate the 
"energy-distance" to each particle; (iii) integrate both the satellite and 
these particles backwards in time in the assumed potential for a Galactic 
lifetime; and (iv) credit the potential with a capture whenever any one of 
these particles is separated from the satellite by a distance dr < 1.8r tide 
and a velocity dv < yj G " r sat ■ 

3. Assign the potential's "score" as the number of successful captures, with 
the most likely potential having the highest score. 
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We tested this algorithm by applying it to 128 stars selected at random 
at the end of each simulation. For each star, 10 test particles were integrated 
in a four-dimensional grid of halo potentials of the form given in equation (JJ) 
with t> circ , q and p varied by 10% around the values in the simulations, and 
x and y axes rotated between 0-90 degrees. The left hand panels of Figure 
|2| show the maximum number of recaptured particles after 10 Gyrs at fixed 
f circ but varying p, q and orientation of the x and y axes, for Models 11 (top), 
13 (middle) and 14 (bottom). The right hand panels show contours of the 
maximum number of recaptured particles in the p — q plane for arbitrary axis 
orientation and v c - ac corresponding to the maximum in the left hand panels. The 
points corresponding to the parameters of the potential in which the simulations 
were originally run are marked with solid squares and the recovered parameters 
are marked with stars. The figure demonstrates that this method is sensitive to 
both the mass distribution and the geometry of the Milky Way. 

We can estimate the error introduced by using N stars to determine the 
potential from the bootstrap method. We create iV(logiV) 2 data sets, each 
of size N, by drawing stars at random (with replacement) from the original 
sample. For each set, we estimate the Milky Way parameters corresponding to 
the grid cell with the maximum score as outlined above. It has been shown (see 
Babu fc Feigelson 1996| and references therein), that the distribution of these 
bootstrapped estimates around the original one will closely follow the intrinsic 
distribution of errors due to sampling. Figure |3] shows the dispersion a w (defined 
as a w = \J (w 2 ) — (w) 2 , where the angles denote averaging) of our bootstrapped 
distribution of estimates for v c i TC , q and p as a function of N. The line in bold 
corresponds to a w oc 1/lOyN. The figure shows that the circular velocity and 
axis ratios of the Milky Way can be recovered to within a few percent using just 
100 stars associated with one of the dwarf spheroidal satellites. 



4. Conclusion: You Can Judge a Galaxy by Its Tail 

In this Letter, we have explored the use of SIM measurements of stars in tidal 
streams as a probe of the Galactic potential. We find that with five-dimensional 
phase-space information for only 100 stars, we can determine the circular 
velocity and shape of the Galactic halo with accuracies of a few percent. This 
measurement would represent more than an order-of-magnitude improvement 
in our knowledge about the Galaxy's mass distribution (Kochanek 1996; Zhao 
1998). 

It should be noted, however, that our discussion has been limited to 



-7- 



constraining four parameters in a specific assumed form for the Galactic 
potential. Ultimately, as with any parameterized inversion algorithm, the 
uncertainty of the method will increase with the number of parameters varied. 
One issue which we have not considered is the degree to which the method 
discussed here can be used to constrain the total extent of the Galactic halo, 
in a manner analogous to that which has been applied recently to tidal tails 
in merging galaxies (e.g. Dubinski et al. 1996, 1998; Springel & White 1998). 
Naively, we expect the measurements proposed here to be most sensitive to the 
mass of the Galaxy enclosed within the region occupied by the debris comprising 
a stream. However, the unusually high precision offered by SIM suggests that 
the influence of the mass distribution on larger scales could be probed if the 
halo is highly flattened. 

The spatial distribution of matter in the Galaxy provides insight into both 
the formation history of our Galaxy and the nature of dark matter. In future 
work we will look at the evolution of tidal debris in lumpy or time-dependent 
potentials and discuss using tidal streams to measure not only the current state 
of the Milky Way, but also to infer its history. 
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Fig. 1. — Top panel — orbital energy relative to the satellite's AE for all particles unbound 
at the end of the simulation of Model 11. Bottom panel — error in the distance estimate 
described in §3.1. The solid lines show the accuracy of distances derived from SIM parallax 
measurements. 
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Fig. 2. — Left hand panels — maximum number of captured particles for fixed t> circ and 
arbitrary q,p and axis orientation as a result of applying our algorithm to Models 11 (top 
panel), 13 (middle panel) and 14 (bottom panel). Right hand panels — maximum number 
of rebound particles contoured in the p — q plane for the most-likely value of t> C i rc identified in 
the left-hand panels. The maxima corresponding to the most likely parameters are marked 
with stars, and the parameters actually used in the simulations are shown with solid squares. 
The contours are spaced by 4 particles. 
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Fig. 3. — Bootstrapped errors in the potential calculated with N stars. The solid triangles, 
squares and pentagons are for Models 11, 13 and 14 respectively. The bold line is given by 
a w = 1/lCh/JV- The dotted line shows the size of one cell of the gridded distribution from 
which a w was calculated. 



